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ABSTRACT 

Based on deep imaging from the Advanced Camera for Surveys aboard the Hubble Space Telescope, 
we present new evidence that stellar feedback created a ~ 1 kpc supergiant HI shell (SGS) and 
triggered star formation (SF) around its rim in the M81 Group dwarf irregular galaxy IC 2574. Using 
photometry of the resolved stars from the HST images, we measure the star formation history of a 
region including the SGS, focusing on the past 500 Myr, and employ the unique properties of blue 
helium burning stars to create a movie of SF in the SGS. We find two significant episodes of SF inside 
the SGS from 200 — 300 Myr and ~ 25 Myr ago. Comparing the timing of the SF events to the 
dynamic age of the SGS and the energetics from the HI and SF, we find compelling evidence that 
stellar feedback is responsible for creating the SGS and triggering secondary SF around its rim. 

Subject headings: galaxies: dwarf — galaxies: individual (IC 2574) — galaxies: irregular — galaxies: 
ISM - stars: formation 

1. INTRODUCTION 

Vigorous (SF) is one mechanism thought to be re- 
sponsible for the creation of the holes and shells ob- 
served in the interstellar medium (ISM) of a wide va- 
riety of galaxies. Following a strong episode of SF, en- 
ergy input into the ISM from young massive stars (e.g., 
stellar winds, Type II supernovae) can provide suffi- 
cient mechanical energy to create these holes and shells 
and may even trigger secondary SF around the edges 
(e.g.. iTenorio-Tagle fc' Bodcnhc imerl Il988t iPuche et al.l 
1992). Although observations of young stellar clusters 
inside these ISM features confirm that this is one possi- 
ble scenario, a majority of holes and shel ls do not con- 
tain bright associated young clusters (e.g.. Ivan der Hulsll 
119961: iRhode et alj|1999t iPasauali et al.ll2008h . One ex- 
planation is that only a select few ISM holes and shells 
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have been observed with sufficient photometric depth 
to detect the associated clusters. Alternatively, a num- 
ber of competing theories suggest the origins of holes 
and shells in the ISM are from mechanisms other than 
stellar feedback e.g., high velocity cloud c ollisions (e.g., 
Kamphuis et al.l [l991h . interactio ns (e.g., iBureau et al.l 
2004h . gamma rav bursts (e.g., lLoeb fe Perna[T l998'l. 
stochastic SF propagation (Ha rris fc Zar itskv 20Q8(). 

Dwarf irregular galaxies (dls) provide excellent envi- 
ronments in which we can study the creation mechanisms 
of these ISM features. Larger galaxies exhibit differential 
rotation and sheer, which can destroy the holes and shells 
in the ISM on timescales as short as 10 7 year. However, 
dls a r e typically solid body rotators (e.g JSkillman et all 
119881 : ISkillmanl[l996l ). which allows features in the ISM 
to grow to large sizes ( ~ 1 kpc) and stay cohe rent over 
longer time scales (e.g.. IWalter fc Brinkslfl999l) . 

Nearby dls are particularly convenient probes of the 
interplay between SF and the ISM, because we can di- 
rectly compare the resolved stellar populations to the 
HI component to better understand mechanisms of feed- 
back . At a distance of ~ 4 Mpc (|Kara chcntsc v et al.l 
l2002f) . IC 2574, an M81 Group dl, contains a "supergiant 
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Fig. 1— (a) HI image l|Walter fc B rinks 1999); (b) optical V-band l|Walter fc Brinks 1993); and (c) color ACS image of the SGS region 
in IC 2574 we use for analysis in this paper. The ellipse in each image is the outline of the SGS and is likely elliptical due to projection 
effects. It has a de-projected diameter of ~ 700 pc. 



shell" (SGS) (deno ted by the re d elli pse in panels (a) 
and (b) of Figure IT] [Walter et al.lll998HWalter fc Brinks! 
1999) and is an excellent candidate for specifically study- 
ing both the formation of the SGS and subsequen t trig- 
gered SF around its rim. IWalter fc Brinks! {1999) clas- 
sify the SGS as an ellipse (~ 1000 x 500 pc), however 
note that many of the holes and shells in IC 2574 are 
elongated due to projection effects (the de-projected di- 
ameter is ~ 700 pc). They further note there is some 
eviden ce that the SGS is a su perposition of two merged 
shells (| Walter fc Brinks! 1 1999f ). which could account for 
its elliptical shape. In addition to the shape, the expan- 
sion o f the SGS is somewhat uncertain. [Walter fc~ Brinks 
(1999) measure an HI expansion velocity of 25 km s 
from a break in the PV-diagram, which is only an in- 
direct measure of expansion, and in fact is consistent 
with stalled expansion as well. We discuss the implica- 
tions of these scenarios in §4. UV o bservations show a 
youn g stellar cluster (~ 11 Myr old, IStewart fc WalteTI 
[20001) inside the SGS, but no Ha emission is observed 
from the interior. Around the rim, strong Ha (and mid- 
IR) emission is obs erved (see panel (c) of Figure [TJ e.g., 
IWalter et ai]|1998t iCannon et ai1l2005f) indicating SF at 
these locations within the past ~ 5 Myr. 

In this paper, we use deep HST/ACS imaging of the 
resolved stellar populations of the SGS to measure its 
recent SFH (< 500 Myr). We use the unique properties 
of resolved blue helium burning stars (BHeBs) to trace 
the spatial locations of the SF episodes, resulting in a 
movie that shows how SF propagated in the vicinity of 
the SGS over the past 500 Myr. From these lines of 
evidence, we are able to directly explore the SF event (s) 
responsible for the creation of the SGS and triggering of 
secondary SF. 

2. OBSERVATIONS AND PHOTOMETRY 

We used HST/ACS to observe the region of IC 2574 
containing the SGS on 2004, February 6 in three wide 
band filters F435W (B), F555W (V), and F814W (I). 
The images were processed using the HST pipeline and 
photometry was extracted u sing DOLPHOT, a version of 
HSTPHOT polphin|[2000l) with an ACS specific mod- 
ule. To account for photometric completeness, we per- 
formed artificial star tests. After filtering out cosmic 
rays, hot pixels, and extended sources from the pho- 
tometry we had 253,736 well measured stars in the ACS 



field of view. Subsequently, we isolated the SGS region 
and found 24,912 well measured stars in the F555W and 
F814W filters. The limiting magnitude of our photom- 
etry in this region is My ~ 0, whic h gives us excellent 
leverage on SF in the past ~ 1 Gyr (jGallart et all 120051 
iWeisz et al.l 120081 ). The photometry was carried out as 
part of a larger study of M81 Group dls and full details 
of the obse rvations, photomet ry, and completeness can 
be found in IWeisz et all (|2008[ ). 

3. THE STAR FORMATION HISTORY 

We measured the SFH of the SGS region (the entirety 
of panel (c) in Figure [ID of IC 2574 from the HST/ACS 
b ased colo r -magn itude diagram (CMD) using the code 
of [Dolphin ( 2002) . This method constructs synthetically 
generated CMDs, u sing the stellar evolution models of 
Marig cTet al.l (|2008h . and compares them to the observed 
CMD using a maximum likelihood technique. To obtain 
this solution, we used a combination of fixed (e.g, binary 
fraction, initial mass function) and searchable (e.g., dis- 
tance, extinction) parameters. For consistency, we set 
all our para meters to the values found for IC 2574 listed 
in Table 2 of IWeisz et all (|2008l ). To quantify the uncer- 
tainties of the SFH, we added the systematic uncertain- 
ties from the isochrones and the statistical uncertainties 
from Monte Carlo tests in quadrature (see IWeisz et al.l 
2008 for details). 

We used the CMD of the SGS region (see Figure [2]) 
as input into the SFH code. The SFH of the SGS re- 
gion over the past 500 Myr is shown in the top panel of 
Figure [31 For reference, we show the SFR of the SGS 
region averaged over the lifetime of the galaxy as the red 
dashed line. We see that over the past 500 Myr, the 
SFR is always higher than the lifetime averaged SFR, 
with notable increases at ~ 200 — 300 Myr and within 
the last ~ 25 Myr. From the most recent time bin (0 — 10 
Myr) we see a dramatic increase in the SFR, which is in 
agreement with the Ha derived SFR for th is region (e.g., 
iMiller fc Ho"dgel 11994 [Cannon et al.ll2005D . Because the 
SGS region, i.e. all of panel (c) in FigureQ] is larger than 
the SGS itself, we tested for the possibility that SF from 
the additional area influenced our SFH. We measured 
the SFH interior to the SGS only and found differences 
which were deemed insignificant. 

3.1. Spatially resolved SFH of the SGS 
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Fig. 2.— The HST/ACS CMD of the SGS region (panel (c) in 
Figure [TJ corrected for foreground extinction. 

BHcBs are powerful tools for studying both the 
ages and locations of SF from ~ 10 — 500 M yr (c.f. 
IDohm-Palmer et al.lfl997l[200llWeisz eTaTEooll . They 
provide us leverage because (a) this evolutionary phase 
is very short relative to the main sequence lifetime 
and (b) BHeBs of different ages do not overlap on 
the CMD (IDohm-Palmer et al.|[l997t iGallart et al.ll2005l; 
iWeisz et al.l l2008n . The result is that we can create 
density maps of BHcBs at different ages, which trace 
the spatial location of SF events over time. Inherently 
uncertainties in distance, extinction, photometry, and 
isochrones can all affect the precise ages of the BHeBs. 
Because of the low metallicities in dls and the high qual- 
ity of our ACS data, we find the observational related 
uncertainties to be small, a nd have calculated isochrone 
based errors to be - 10% (jWeisz et al.ll2008t ). In addi- 
tion to uncertainties in ages, the locations of the BHeBs 
can be influenced by proper motions of the stars. The 
stars we consider are thought to belong to clusters and 
are thus bound by gravity. Although they will not re- 
main gravitationally bound indefinitely, the time scale 
for dissip ation of such stellar associa tions is ~ 600 Myr 
to 1 Gyr (|Dohm-Palmer et al.| [T997), much longer than 
the time scale we consider for the SGS. To transform the 
density of BHcBs to a SFR per area, we multiply each 
density map by the corresp onding SFR, measured using 
the code of iDolphinl (|2002f) . divide by the image area, 
and interpolate in both space and time to create a spa- 
tially resolved SFH movie. The final spatial resolution of 
the movie is ~ 8", similar to that of the HI observations, 
and the movie frames are then interpolated onto a grid 
with 5 Myr sampling. 

Figure [4] shows selected still frames from the spatially 
resolved SFH movie of the SGS region; the outline of the 
SGS is indicated by the white ellipse. We can see that 
in the oldest frames, 500 and 350 Myr, there is no major 
SF going on in the region. From 200 — 300 Myr, we 
see that there is an elevated period of SF in the center 



Fig. 3.— (top) The SFH of the SGS region over the past 500 Myr. 
The red dashed line indicates the SFR of the SGS averaged over 
the lifetime of the galaxy; (bottom) The cumulative energy from SF 
in the SGS over the past 500 Myr calculated with STARBURST99. 

of the SGS. A period of relative quiescence lasts until ~ 
100 Myr, when a SF event begins to strengthen inside 
the central regions of the SGS. The peak strength of this 
central event is at ~ 25 Myr, after which SF subsides 
centrally, but increases in intensity on the rim. 

The location of the SF events in the center of the SGS 
and on the rim are in excellent agreement with obser- 
vations from other wavelengths. Th e central event co- 
incide s well with UV observations by I Stewart fc Walter I 
(2000). The triggered SF o n the rim is traced by Ha 
and Spitzer IR imaging (e.g.. ICannon et aill2005D and is 
in excellent agreement with our 10 Myr frame. The 200 
— 300 Myr event is beyond the look back time of the 
UV and IR fluxes, and is only traceable with the deep 
photometry presented here. 

4. DISCUSSION 

To understand the current state of the SGS, we con- 
sider the energetics of both the HI and SF, the role of 
triggered SF, and the timing of the events that led to the 
evacuation of the HI and secondary SF around the rim. 
To do this we examine the efficiency of energy, i.e. the 
feedback efficiency, transfered from stellar feedback into 
moving the ISM in the context of both an expanding and 
stalled shell. 

Following the model of lChevalierl (|1974f ). one can com- 
pute the energy necessary to evacuate the HI mass inte- 
rior to the SGS to be ~ 2x l 53 erg (e.g.. IWalter et all 
fl998t IWalter fc Brinks! fl99l ICannon et all l2005h . Us- 



ing STARURST99 (|Leitherer et alJll999T T"we can inde- 
pendently calculate the energy associated with SF in the 
SGS region by inputting our SFH and assum ing the same 
IMF as was used to measure the SFH (see IWeisz et al.l 
2008). We computed the total cumulative kinetic en- 
ergy (i.e., stellar winds and SNe) of the SGS region over 
the past 500 Myr to be ~ 3xl0 55 erg (Figure ©. Be- 
cause of our relatively coarse time resolution, we need to 
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Fig. 4. — Selected still frames from the spatially resolved recent SFH of the SGS region. The white ellipse corresponds to the elliptical 
outline of the SGS itself shown in Figure Q] The spatial resolution of the images is ~ 8", similar to that of the HI observations. The movie 
can be seen at http://www.astro.umn.edu/~dweisz/2574/sgs.mov . 



consider that our elevated SF periods, may be shorter 
bursts averaged out over a longer time. Thus, we used 
STARBURST99 to model the SFH as both a series of 
constant episodes of SF, as well instantaneous bursts for 
the elevated SF episodes and found the energetics to be 
consistent. 

At first pass, it may seem that SF generates an un- 
reasonably large amount of energy compared to the HI 
energy requirement. However, by considering the feed- 
back efficiency we can reconcile the apparent discrepancy. 
Assuming that each SN imparts 10 51 erg in into mov- 
ing the ISM we find a ~ 1% stellar feedback efficiency, 
which is on the low end of the spectrum that ranges 
from ~ 1 - 20% (e.g-lTheis et al.lll992t iCole et aLl ll99l 
iPadoan et aDl!997t rThornton et alJll998h . However, we 
have really computed a lower limit on the efficiency, as 
the HI energy is a lower bound, and the SF energy is 
an upper bound. The HI energy was computed using 
an assumed smooth uniform HI distribution inside the 
SGS. Given that the central SF event in the SGS was ex- 
tremely powerful, it is likely that there was a higher den- 
sity o f gas interior to the SGS in th e past (e.g . .lOey et al.l 
l2002t) . For example, the model of I Chevalier] (|1974h pre- 
dicts that an increase in a factor of 10 in the gas density 
would change our feedback efficiency to ~ 10%. When we 
consider the energy input into the ISM due to SF, only 
a fraction of this energy is used to move the ISM, with 
the rest lost to heating and radition. Simulations sug- 
gest ~ 10 50 erg per SN is a more realist ic amount of SN 
energy imparted into moving the ISM (|Cole et al J 11991 
iThornton et al.|[l998f ). which also increases our feedback 
efficiency. Finally, if the SGS has stalled expansion and 
broken out of the plane of the galaxy, the energy from SF 
simply escapes and does not go into moving the ISM. De- 
spite the large uncertainties in the calculations of both 
energy quantities, we find that our feedback efficiency 
falls in the range of the values typically calculated from 
simulations, indicating the SF generated an appropriate 
amount of energy for creation of the SGS. 

The SGS presents a clear case of not only sequen- 
tial SF, but also triggered (i.e., causally connected) SF. 
From the spatially resolved SFH movie (Figure 2]) , we 



see clearly that the central burst predates the SF on the 
rim, which is in excell ent agreement with age estimates 
in the literature (e.g.. IWalter et ail 119981 : ICannon et al.l 
[20051: IPasauali et at] 12008). Establishing a causal rela- 
tionship between SF events is generally non-trivial and 
only a han dful of cases of tr iggered SF have been estab- 
lished (see lOev et al"1 120051 and references therein) . In 
the case of the SGS, we have a recent SF episode interior 
to the SGS, which peaked ~ 25 Myr ago and the bulk 
of the recent SF on the rim is < 10 Myr old, indicated 
by both the SF movie and the Ha observations. Thus, 
it appears that the central SF event initiated the expan- 
sion of the shell, which swept up the ISM, and when the 
gas on the rim becomes sufficiently dense, secondary SF 
began. As a slight twist to this scenario, we can specu- 
late about the role of the older (200 — 300 Myr), central 
SF event. Because of the location of the older event and 
energy input into the ISM, it is possible that this older 
event cleared out some or all of the SGS, and the more 
recent event simply served to accelerate the expansion 
and trigger the SF around the rim. Or perhaps the older 
event created one spherical shell, which merged with a 
younger shell to account for the ellipticity of the SGS. 
To better our understanding of the role of both young 
and old SF events in reshaping the ISM, we will need to 
employ similar analysis techniques on a statistically sig- 
nificant sample of holes and shells in the ISM of nearby 
dls. 
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